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ABSTRACT: 

An investigation of elastic buckling of gun-launched motor cases 
is considered. The SATANS finite difference computer program was 
used for the analysis. Comparative results from other buckling 
models were obtained to verify the results. The results are sum- 
marized in a set of design curves and show that for the proposed 
design parameters, the elastic behavior is exceeded before buckling 
can occur. 
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NOTATION 



a, a m acceleration, maximum acceleration (in g's) 

E Young's modulus of elasticity 

G Shear modulus 

h thickness of motor case 

L length of motor case 

N buckling load per unit circumference 

p.j radial pressure due to slumping propellant 

p Q escaping breech force pressure 

P s motor case body force due to acceleration 

P a aft edge force 

Pf forward edge force 

R radius of motor case 

Wf weight of structure fore to motor case 

v Poisson's ratio 

a axial stress 

Oy yield stress 

<d density 
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INTRODUCTION 



This report is concerned with the elastic buckling of gun- launched 
motor case structures which are subjected to large axial g loading 
during launch. Acceleration loading up to 10,000 g is typical. The 
axial loading is accompanied by a lateral loading due to an internal 
slumping propellant which also provides additional stiffness. In a 
previous report (ref. 1, Appendix D) a computation of the buckling load 
for a motor case with 0.08 inch thickness, 1.5 inch radius, and 9.5 inch 
length showed buckling occurs at a stress of 356 ksi , exceeding the 
elastic limit of the high performance heat-treated steel (200-250 ksi). 
Although that analysis indicated elastic buckling was not probable, a 
parametric study of the elastic buckling of the case for various values 
of radius, length, thickness and forward mass was undertaken in the 
event steels with higher heat treatment were used. A finite difference 
computer program, known as SATANS (ref. 2) and capable of performing 
geometrically nonlinear static and dynamic analyses of arbitrarily 
loaded shells of revolution, was used to compute the static buckling 
loads. 
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DESCRIPTION OF THE MOTOR CASE 



The motor case structure is a cylindrical shell of 4130 heat- 
treated steel alloy. From MIL-HDBK-5A, February 8, 1966: 

6 

Elastic modulus in tension and compression, E = 30x10 psi 

6 

Modulus of rigidity, G = 11x10 psi 
Poisson's ratio v = 0.32 
Density, u> = 0.283#/in 

The design parameters furnished by The Naval Weapons Center, 
China Lake, are: 

radius, R = 2.5 inches 

length, L, between 19 and 24 inches 

thickness, h, between 0.125 and 0.250 inches 

forward structure weight, W^, between 45 and 70 lbs. 
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DESCRIPTION OF THE EXTERNAL LOAD 



The external load on the motor case structure was originally 
idealized as shown in the accompanying figure. 

is the force per unit circum- 
ference along the forward edge 
and is associated with the 
acceleration of the forward 
structure. 

A 

P f = 2 ttR * a m 

where a m is the maximum 
m 

acceleration in g's. 

P s is the axfal distributed force 
per unit circumference due to 
the acceleration of the motor case itself, 

p = coha 

s m 

where go is the case material density. 

P, is the aft force per unit circumference obtained by equilibrium 

a 

considerations 

P a ■ P f + p s ,L 

Finally, p. and p Q are the lateral pressures associated with the 
slumping propellant and escaping breech pressure, respectively. These 
two pressures were neglected in the analysis. The breech pressure p Q 
can be minimized by good obdurator design, and it would be difficult 
to determine the magnitude under any circumstances. The effect of the 
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slumping propellant pressure p. was examined in detail and is the sub- 
ject of another report. It was determined that the maximum internal 
pressure would be of the order of 10 ksi for a motor case with a 2.5 
inch radius, a 10,000 g acceleration, and a 20 inch length. Theory 
predicts and preliminary analyses verified that the internal pressure 
stiffens the cylinder. This fact, together with the uncertainty of 
the propellant parameters, provided the basis for neglecting the 
internal pressure in the analyses. 

ANALYSIS 

The validity of the results, which are given in the next section, is 
considered here. The boundary conditions on both ends are taken as 
simple supports. 

The various thin shell theories are associated most often with 
slightly different approximations and assumptions. For example. Love, 
Donnell, Reissner, Flugge, Sanders and others have thin shell theories 
associated with the particular thin shell models they have proposed. 

The finite difference program SATANS used in this analysis is based 
on the Sanders model for thin shells (ref. 2), which includes the non- 
linear terms in the strain-displacement and equilibrium relations. The 
results obtained from SATANS for a cylinder subjected to an axial load 
are compared to those from three other sources; the Flugge equations 
(ref. 3), the classical (linear) buckling equation (ref. 4), and 
empirical results (ref. 5). The results for the critical buckling 

load N and corresponding stress o are given in Table 1 for various values 

6 

of R/h. In all cases. Young's modulus and Poisson's ratio are 30x10 psi 
and 0.32, respectively. 
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Table 1. 



N is in lb. /in. 



o is in psi . 





SATANS 


FLUGGE 


CLASSICAL 


EMPIRICAL 


R/h = 288 


N * 360 


N = 315 


N = 439 


N = 144 


L = 40", R = 2" 


o = 52,700 


a = 45,500 


a = 64,300 


o = 20,800 


R/h = 144 


N = 1062 


N = 1112 


N = 1755 


N = 633 


L = 40", R = 2" 


o = 76,600 


a = 80,400 


a = 126,800 


a = 45,700 


R/h = 92 


N = 3264 


N = 3000 


N = 4386 


N = 1656 


L = 40", R = 2" 


o = 149,000 


a = 137,000 


a = 200,000 


o = 75,600 


R/h = 16 


N = 46,300 


N = 43,300 


Not 

Applicable 


N* = 58,500 


L = 40", R = 2" 


a = 370,000 


a = 346,000 




o* = 468,000 



In all cases the results given in Table 1 show the classical (linear) 
buckling theory gives buckling loads which are not in agreement with 
any of the other sources. Its inclusion here was simply to demon- 
strate this shortcoming. It is also noted that for R/h = 16 the 
classical theory cannot be used. 

The first three cases of the table show that for very thin shells, say 
R/h > 50, the SATANS and Flugge theories agree quite well but give 
buckling loads twice as large as the empirical results. This differ- 
ence has been noted by many investigators (see ref. 5, page 509), and is 
accounted for by shell imperfections and load eccentricity. For the motor 
case dimensions provided by NWC, the R/h ratios lie between 10 and 20. 
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Here Flugge and SATANS agree quite well and if the empirical results are 
extrapolated, additional agreement is obtained. In addition to this 
corroboration between theory and empirical results, it is noted that for 
these shells with R/h = 15, the theoretical solution is on the safe side 
of the extrapolated empirical result. The superscript asterisks on the 
empirical N and o for R/h = 16 indicates these results have not been 
achieved in the laboratory, rather they have been obtained by extrapo- 
lation to provide some comparison here. 

DESIGN BUCKLING CURVES 

The results obtained from SATANS for the range of design parameters 
given by NWC are shown in figures 1 through 5. Figures 1 through 3 are 
plots of acceleration versus L/R for the minimum and maximum h/R ratios 
(0.05 and 0.1 respectively). Each figure is for a specified forward 
weight, W^, of 45, 60, and 75 pounds respectively. The buckling curves 
are specifically labeled. The remaining curves denote failure by yield- 
ing for 200, 250, and 300 ksi steels. The most Important observation 
here is the one anticipated by the previous analysis; elastic buckling 
will not occur. If buckling is to occur, it must be plastic. 

The small variation of acceleration with L/R ratio suggests the alternate 
plots of figures 4 and 5 which show acceleration versus for minimum 
and maximum h/R ratios for fixed L/R ratios. 

DESIGN TECHNIQUE 

Figures 1 through 3 or figures 4 and 5 may be used for design by select- 
ing any three out of the four quantities, a, h/R, L/R, and W f , and 
obtaining the remaining fourth quantity from the lower of the two curves 
associated with buckling and elastic yield. In all cases considered 
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here the latter is the failure mode. 

The values of h/R equal to 0.05 and 0.10 are the minimum and maximum 
values of h/R obtained from the NWC design parameters, R = 2.5 inches and 
h between 0.125 and 0.250 inches. The buckling curve for h/R equal to 
0.075 was determined to provide a means of interpolation for values of 
h/R intermediate to 0.05 and 0.10. The proximity of the 0.075 curve to 
the mean of the 0.05 and 0.10 curves suggests that linear interpolation 
may be used. As an example consider the case for h/R = 0.625, L/R = 7.6 
and Wf = 60 pounds. From figure 4 we find a m equal to 17,600 g and 
42,300 g for h/R equal to 0.05 and 0.10 respectively. 3y linear inter- 
polation the maximum acceleration for static buckling when h/R = 0.0625, 
L/R = 7.6, and Wf = 60 pounds is 

a m = 6 /5 I *|o (17,600 + 42,800) = 35,200 g. 

Because figures 1 and 5 have log scale ordinates, direct scaling of the 
curves is not valid. 

The relation governing the initial yield curves is 

F ma ( W f + W s>* 
a ~ A ~ A " 2*Rh 

where W s is the weight of the motor case, Wf, is the weight of the fore 
structure, and the acceleration a is in g units. Solving for the accel- 
eration when the stress is at yield gives 

_ 2irRh _ 2nRh 

a m “ Wf + W s a y Wf + 2nRh . <d ¥ 

This equation may be solved for any values of the design parameters. 

It may be noted that among all variables the maximum acceleration, 

a m , is linearly proportional to yield stress, o y , only. 
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